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We image the flow of a nearly random close packed, hard-sphere colloidal suspension (a 'paste') in 
a square capillary using confocal microscopy. The flow consists of a 'plug' in the center while shear 
occurs localized adjacent to the channel walls, reminiscent of yield-stress fluid behavior. However, 
the observed scaling of the velocity profiles with the flow rate strongly contrasts yield-stress fluid 
predictions. Instead, the velocity profiles can be captured by a theory of stress fluctuations originally 
developed for chute flow of dry granular media. We verified this behavior both for smooth and rough 
boundary conditions. 
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Understanding the deformation and flow, or rheology, 
of complex fluids in terms of their constituents (colloidal 
particles, polymer chains or surfactant aggregates) poses 
deep challenges to fundamental physics, and has wide 
industrial applications T. The experimental study of 
complex fluid rheology typically starts in a rheometer, in 
which stresses and strains are applied and measured in 
well-defined, 'rheometric' geometries ('cone-plate', etc.). 
Translating rheometer data to more complex flows is non- 
trivial, but well developed in polymers (see, e.g., |2|). 

The understanding of colloidal flows lags consider- 
ably behind and despite their equal practical impor- 
tance 11] , studies on model systems have been carried 
out only recently [5]. Compared to polymers, colloids 
pose some unique challenges. Concentrated suspensions 
('pastes') are generally non-ergodic (or 'glassy'), so that 
any flow involves non-linearities (e.g. yielding [6l [7] or 
shear thickening [Hlin])- Moreover, specific geometries in 
applications may involve dimensions comparable to sin- 
gle particles and lead to confinement effects, such as in 
micro-fluidics [10]. The most quantitative theory for qui- 
escent colloidal glasses, mode coupling theory, has only 
recently been extended to deal with simple shear [11 . 

Here we present an experimental study of the flow of a 
hard-sphere suspension at nearly random close packing, a 
'paste', in a twenty-particle- wide square capillary. Pastes 
are ubiquitous in industry, where their unique rheology 
presents many challenges and opportunities [3] . The sim- 
plicity of the geometry is appealing from the fundamental 
perspective. It also makes direct contact with microflu- 
idic applications [10] . Using fast confocal microscopy, we 
tracked the motion of individual colloids and measured 
the velocity profiles in channels with both smooth and 
rough walls. Despite the colloidal nature of our suspen- 
sion, we find significant similarities with granular flow, 
itself of wide applied [T^ and fundamental fTT interest. 

We used sterically stabilised polymethylmethacrylate 
(PMMA) spheres of diameter D 2.6 ± 0.1 ^m 
(from confocal microscopy) fluorescently labelled with ni- 
trobenzoxadiazole, and suspended in a mixture of cyclo- 
heptylbromide and mixed decalin (viscosity 2.6 mPa-s) 
for buoyancy matching at room temperature T = T^. 



Their Brownian time is tb — /24:Dq ~ 4.4 s with Dq 
the bare diffusion coefficient. 
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FIG. 1: Velocity profiles in the lower half of the capillary 
(width 2a = 50 jj,m) with smooth walls reconstructed in 3 jj,m 
steps. The average velocity (V) is in units of the average 
velocity of the central, unsheared, 'plug', (Vc) ~ 20 /ims"^. 

A dense suspension (volume fraction > 0.63, from 
confocal microscopy) was obtained by centrifuging (at 
T > Tr to remove buoyancy matching). A constant 
pressure gradient, Vp, is applied to drive the suspen- 
sion into a square borosilicate glass micro-channel (Vit- 
rocom Ltd; side 2a — 50 /im) [TJ], whose inner walls 
were either untreated and smooth, or coated with a dis- 
ordered monolayer of colloids and thus rough on the par- 
ticle level. The coating particles, which are sightly larger 
(-Dcoat = 2.8 /xm) and more polydisperse PMMA spheres 
(right inset to Fig. |3]), were applied by flUing with a 
dilute suspension and attached by heating in a vacuum 
oven (110 °C). 

The flow (along x) across the full width of the chan- 
nel {\y\ < a) was imaged with a Visitech VTeye confocal 
scanner in a Nikkon TE Eclipse 300 inverted microscope. 
We collected 44 //mx58/im images (107 frames per sec- 
ond) at depths —a < z < +0.2a (z = —a is the lower sur- 
face). From the two-dimensional images, we located par- 
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FIG. 2: Averaged velocity profile as a function of 2y/aefi for 
smooth (o) and rough (V) walls at z = 17 /xm. Full lines: fits 
from integrating Eq. [7] dotted lines: extents of sheared zones. 
(Vc) is 114 /ims~^ for the smooth wall case and 80 /ims~^ for 
the rough wall case. Inset: width of shear zone in units of 
particle diameter as a function of AV. A and T: two runs for 
rough walls; • and ■: two runs for smooth walls. 
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FIG. 3: Normalized velocity profiles at 2; = 17 /im versus 
2j//acft . We use the profiles fitted from integrating Eq. [t] the 
error bars show the spread in the measured data. Left inset: 
normalized transverse fiuctuations of individual particles in 
the shear zone for the two boundary conditions, versus nor- 
malized time I/tb or tocoZ accumulated strain t-d{V{y)} /dy = 
t ■ 7(y), top axis. The traces show larger fluctuations in the 
rough wall case. Right inset: confocal microscopy snapshot 
of the flow at 2 = 17 fim for a channel with rough walls; the 
particles marked with • are attached to the walls. 



tides with a resolution Sx,5y ~ 50 nm [15J. A key step 
for correct tracking is to follow colloids in a 'comoving' 
frame [14j . This involves removing from the 'raw' coordi- 
nates the advective motion Ax{y), which is obtained as 
the a;-shift that maximizes the cross correlation between 
image strips around y in successive frames. From time- 
dependent coordinates (restored in the lab frame) , we ob- 
tained particle velocities |25j . After start-up transients, 
we found oscillations in the particle velocities (data not 
shown) at 0.1 Hz (slow flow) to 1 Hz (fast flow). This 
feature, which we discuss elsewhere, is ubiquitous in the 
pipe flow of pastes [5l[l6]. Here, we restrict ourselves to 
steady-state velocity profiles, {V), obtained by averaging 
over 20 s (slow flow) to 10 s (fast flow), corresponding 
to ~ 2000 to ~ 1000 frames, respectively. We imaged 
at a; ~ 0.5 cm from the entrance to the capillary (corre- 
sponding to ~ 2000 particles), where entry effects have 
died out and the results show negligible x dependence. 

Typical data for smooth walls are shown in Fig. [T| 
For z > 10 /im, each velocity profile consists of a shear 
zone close to the walls and a nearly unsheared central 
'plug'. This plug shrinks at smaller z, i.e. closer to 
the bottom wall. We also observe wall slip with ve- 
locity (Vs). The profiles for smooth and rough walls. 
Fig. [2] are qualitatively similar, but the latter displays 
considerably smaller wall slip and larger shear zones. 
Note that for rough walls, we use an effective half width 
aes = a — Dcoat — D/2; for smooth walls, UcS = a — D /2. 

The dependence of the velocity profiles on the overall 
flow rate is striking. We define the width of the shear 
zone h as the distance from the wall where the flow speed 



is (14) -I-0.95AF where = (K) - (V;) is the difference 
between the averaged center and wall speeds (with AF/a 
the average shear rate). We see in the inset to Fig.[2]that 
h is independent of AT^, i.e. the central plug remains 
essentially constant for flow rates varying by more than 
one decade. Associated with this, the normalized velocity 
profiles {{V) — {Vs))/AV for different flow rates collapse 
onto two different master curves for rough and smooth 
walls. Fig. [3] 

Let us now compare the observed behavior to what 
is expected from bulk rheology. Considering the system 
as a glassy yield stress fluid, the stress r versus strain 
rate 7 relation typically displays Herschel-Bulkley shear 
thinning behavior: r — Tyidd oc 7" (n « 0.5) [6] with 

'^yield 

the yield stress. Plug flow is ubiquitous in all yield stress 
fluids, and has been extensively studied [T7]. Thus, our 
observation of plug flow is, in itself, unsurprising. In yield 
stress fluids though, if the maximum applied shear stress 
in the channel exceeds Tyidd, we have plug flow with shear 
zones whose boundaries are given by r = Tyieid; hence, 
the width of the shear zone, b, grows with the shear rate 
and eventually leaves a vanishingly small 'plug' '17]. This 
is not what we observe. 

A clue to what may be happening in our experiments 
comes from the observation that at the typical strain 
rate encountered in the shear zones, Tg^AV/a) ^ 50, 
conventional rheology would lead us to expect severe 
shear thickening. Shear thickening in pastes is very 
far from completely understood on the microscopic level 
[5J UHl Hn]; but we may speculate that for colloids sta- 
bilised by short grafted polymers, interparticle friction 
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may become important as they jam against each other 
driven by shear. Indeed, the measured transverse fluc- 
tuations of particles clearly show the motions induced 
by them 'bumping' along neighboring layers, Fig. 3 in- 
set. We are thus led to consider analogies with friction- 
dominated granular flow. 

Indeed, both the width of the shear zone b ^ 6D 5D 
for smooth walls) and the lack of flow rate dependence, 
are strongly reminiscent of observations in gravity-driven 
'chute flow' of dry granular materials |20j . In particular, 
Pouliquen and Gutfraind observed b « 6D in the two- 
dimensional flow of discs (diameter D) for channels of 
widths lOD < 2a < 44Z? and developed a model to ac- 
count for their observations [^Hl [H]. Below we extend 
their model to three dimensions, and show that it indeed 
predicts b independent of AV. 

Following [20] we start with the components of the 
stress tensor r in the fully-developed flow of a continuum 
medium in the x direction along a pipe, which satisfies 
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where tq = SaVp/n'^. This constant sets the scale for the 
(z-dependent) maximum stress at the wall, r,nax- 
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FIG. 4: (Top) Fitted fluctuating stress amplitude At in units 
of Tniax as a functiou of the 'net shear', AV; symbols as in 
Fig. [2] (Bottom) Schematic of the fluctuating shear stress su- 
perimposed on the continuum T^cy The dashed line shows the 
quadratic fit used for the integration. In the zone where the 
fluctuating stress overcomes Tyieid, flow results. 



To determine the normal stress, Tyy we impose a 
Coulomb friction condition at the walls: 



Tyy — Aiwall'''max j 



(4) 



where /iwaii is the friction coefficient between the sus- 
pended particles and the wall; this constant Tyy satisfies 
Eq.[2][26]. Following the Coulomb criterion, the material 
yields if r^y > Tyidd = MbuikTyy, with //bulk the friction 
coefficient inside the material. We expect /iwaii £ A*buik, 
with '=' for a rough wall coated with particles, and '<' 
for smooth, glass walls. From Eq. |4]we have: 



T'yicld — (A^bulk/MwalOT'n 



> T„ 



(5) 



Since r^y < Tyidd, the bulk never yields and the whole 
material slips as a plug. Note from Eq. [5] that ryjeid 
increases with Tmax (and therefore AV). This reflects 
the rising normal stress Tyy, Eq. |4] which increases the 
friction between particles and thus TyidH [27]. 

We now assume the presence of stress fluctuations 
which, when added to the continuum T^y, Eq. [3] may 
take the local stress above Tyidd, Fig. H Following [20] 
we use a 'Boltzmann' ansatz for the yielding probability: 



Pyidd oc exp 



Tyieid (^)- I T^yjy.z) I 

At 



(6) 



where At = At(z) is the amplitude of stress fluctua- 
tions, taken to be independent of y at any particular z. 
The model is completed by the simplest possible ansatz 
relating shear rate to Pyidd: dV{y,z)/8y cx Pyidd [28| . 
At a fixed z = zq, therefore, we have 



dV{y,zo) 
By 



constant x exp 



ixy 



Ai 



(7) 



In order to perform the integration analytically we fitted 
the stress profile, Eq. [3| with a quadratic approximation; 
Fig. [4] (bottom) shows both the analytic expression and 
the fit. We then substituted the fitted stress into Eq. [7] 
integrated and fitted the resulting velocity profile to the 
measured data. Example results are shown in Fig. [2] 

The fitted stress fluctuation amplitudes normalized to 
the wall stress, AT/Tmax, at z = 17 /im for both smooth 
and rough walls are plotted against AV in Fig. |4] As 
7'max increases with flow rate. At increases proportion- 
ally so that their ratio remains constant. Combined with 
Eqs.|3]and Sjthis means that {Tyieid — \'''xy\)/AT in Eq. [6] 
is independent of Tmax, resulting in a constant value for 
6, Fig. |2] and collapse of the velocity profiles. Fig. [3] 

The fact that At scales as the wall stress T^ax suggests 
that the stress fluctuations are controlled by what hap- 
pens at the boundaries. Two dimensional simulations 
of dry granular matter indeed relate stress fluctuations 
to inhomogeneities in the friction-dominated force chains 
resulting from contacts at the walls [23]. Here, we find 
(for z ^ a) that AT/Tmax for rough walls is somewhat 
larger than for smooth walls, Fig.|4| this difference is di- 
rectly correlated with the observed difference in 6, Fig. [2] 
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The presence of such a difference, and its sign, is unsur- 
prising. Flow along a smooth wall is less 'bumpy' and 
can thus be expected to generate a lower level of stress 
fluctuations which propagate through the system, as also 
shown by the two traces in the left inset to Fig. [s] ^ . 
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FIG. 5: Average density profile aX z = VT fim for rough walls, 
in units of the density ric in the center. Dashed lines: average 
width b of the shear zones in channels with rough walls. 

The analogy with granular flow extends further than 
the constant shear zone width. In [20j , a variation of 
the density across the channel was observed, the 'plug' 
being ^ 10% denser than the edge of the shear zone. We 
observe the same feature in our channel flows. Figure [5] 
shows a typical density profile ai z — 17 /im, measured by 
counting particles per unit area and averaging 10* frames; 
this profile is essentially independent of the overall flow 
rate. The observed density reduction in the shear zone 
is not unexpected, since dilatancy is required in flowing 



particulates (dry or wet) at such high packing fractions. 

Our findings also have implications concerning the role 
played by confinement. Since D/2a ~ 20, one might 
expect strong size-dependence. However, the width of 
the shear zone b we observed, is always smaller than the 
width of the channel, a and size effects should only be- 
come significant for channels with a ^ b. 

To summarize, we have measured the flow properties of 
a nearly-random-close-packed hard-sphere colloid driven 
by a constant pressure gradient in a twenty-particle- 
diameter square channel by tracking the motion of in- 
dividual particles using fast confocal microscopy. At all 
flow rates, we observed a central, almost unsheared plug, 
and peripheral shear zones. In contrast to the prediction 
of yield-stress fluid rheology, the size of the shear zone 
remained constant as the flow rate was increased. We 
explained this by appealing to a model originally set up 
for the gravity-driven 'chute flow' of dry granular mate- 
rials [50]. This model should be applicable in our case 
if the yield stress, Tyidd, is dominated by inter-particle 
friction. The model predicts that stress fluctuations can 
bring about yielding even when the average stress is be- 
low Tyioid- Quantitative fits to our data were obtained. 
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